Excitons in long molecular chains near the reflecting interface 
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We discuss coherent exciton-polariton states in long molecular chains that are formed due to the 
interaction of molecular excitations with both vacuum photons and surface excitations of the neigh- 
boring reflecting substrate. The resonance coupling with surface plasmons (or surface polaritons) of 
the substrate can substantially contribute to the retarded intermolecular interactions leading to an 
efficient channel of the decay of one-dimensional excitons with small momenta via emission of surface 
excitations. The interface also modifies the radiative decay of excitons into vacuum photons. In an 
idealized system, excitons with higher momenta would not emit photons nor surface waves. For a 
dissipative substrate, additional exciton quenching takes place owing to Joule losses as the electric 
field of the exciton polarization penetrates the substrate. We discuss how these effects depend on 
the polarization of molecular excitations, their frequency and on the distance of the chain from the 
substrate. 

PACS numbers: 78.67.-n, 71.35.Aa, 71.36,+c, 73.20.Mf 
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I. INTRODUCTION. 

The interaction of the electromagnetic field with 
a molecular excitation in an aggregate of identical 
molecules leads, on the one hand, to the delocalization of 
the excitation over the aggregate and, on the other hand, 
to a modification of its radiative decayi Both notions of 
excitons and polaritons are used in the literature on such 
delocalized excitations. In crystalline structures, the ex- 
citation can be spatially coherent and then it is character- 
ized by its wave vector as a proper quantum number. In 
this paper we discuss how electric-dipole-active coherent 
excitations in linear crystals are affected by the presence 
of a neighboring metallic/dielectric half-space. 

A variety of one-dimensional (ID) electronic systems 
available nowadays, such as conjugated polymers, J— 
and H— aggregates, semiconducting quantum wires and 
carbon nanotubes, exhibit interesting optical properties 
and are considered for potential applications in optoelec- 
tronics; their spectroscopy is an active research area. 
Successes at the synthesis and fabrication of these sys- 
tems have resulted in the continuously improving quality 
and the increase of their "conjugation length" L which 
may exceed the appropriate electromagnetic wave length 
A. Perhaps one of the most noteworthy achievements in 
this regard is a recent observation^ of a macroscopic co- 
herence of a single exciton in polydiacetylene chains of 
L ~ 10 fj,m that allowed to discuss an issue of "an ideal 
ID quantum wire" It is also relevant to note a physically 
very different but conceptually related class of excitations 
in chains of "dipole-coupled" nanoparticles (see, e. g., a 
recent Ref. and multiple references therein) studied for 
photonic and plasmonic applications. 

Long before the modern experimental advances, it was 
shown— £ that the coherent interaction of low-dimensional 
(ID and 2D) excitons with the electromagnetic field is 
drastically different than in 3D systems. The radiative 
decay of low-dimensional excitons is strongly enhanced in 



the region of their wave vectors |q| < k, where k = uj/c 
(c is the speed of light in vacuum) is determined self- 
consistently by the excitation energy E(q) — Tilo, while 
the excitations with q| > k would not radiate, as is re- 
quired by the energy and momentum conservation for 
an exciton-photon system. More specifically, for ID ex- 
citons in the molecular chain in vacuum, the radiative 
width r = K/t (r being the decay time) depends on the 
wave vector q as 
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where p is the magnitude of the molecular dipole transi- 
tion moment p that makes angle 9 with the chain axis, 
a the intermolecular spacing along the chain, and the 
step-function. As compared with the radiative width 



r°(fc) = Ap 2 k 3 /3 



(2) 



of a single molecule, Eq. {1} exhibits an enhancement fac- 
tor of ~ A/a = 2-k jka ~ 10 2 — 10 3 in the optical region of 
the spectrum ( "superradiant" states) . This enhancement 
has been discussed in the context of various systems (see, 
e. g., Refs. 601111 and references therein). Equation (p} 
also shows how the polarization of the transition dipole 
affects the q— dependence of the decay rate. 

Importantly for applications, it is possible to manip- 
ulate the optical properties of molecular excitations and 
to form new hybrid excitations by putting molecules or 
molecular aggregates in the vicinity of interfaces and in 
dielectric microcavities^ Well-known, e.g., is an oscil- 
lating dependence of the radiative width on the distance 
of a single molecule from the planar interface resulting 
from the interference of the radiative fields of a molecular 
transition and its image i 10 ' 11 The dipole-dipole interac- 
tion gets also modified in the vicinity of a surface or in 
the cavity i 12 i 13 ' 14 In the case of a linear molecular crys- 
tal, the environment can lead to qualitatively interesting 
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FIG. 1: Decay width F(k,q) of ID exciton-polaritons as a function of the reduced wave vector q/k for a fixed value of 
k — 27r/400 a and substrate's dielectric constant e = —10. Two rows correspond to different dipole-to-image-dipole distances 
d: upper panels (a-c) are for d = 600 a, lower panels (d-f ) for d = 10 a. All distances are measured in units of intermolecular 
spacing in the chain, a. Each of the trio of the panels in a row corresponds to different exciton polarization, indicated by the 
boldface letters. The chain is parallel to the x-axis and situated above the substrate whose surface is the xy-plane. Decay 
width is shown with respect to T v = 2ivp 2 k 2 /a which is the radiative width of the x-polarized exciton in vacuum at q = 0, 
Eq. {T}. The overall vacuum benchmark results {TJ are shown with the dash-dotted lines. The total decay width is displayed 
by thick solid lines, the dashed solid lines (when distinguishable) show the part of the width due to the decay into vacuum 
photons only. The vertical dash lines indicate the position of the surface plasmon wave vector, K p /k, Eq. Q. 



coherent effects as it is the interaction of many molecu- 
lar transition dipoles (and their images) that would de- 
termine the properties of the excitation. Philpott,— for 
instance, pointed out that by placing a linear chain near 
a transparent substrate, one could probe some of oth- 
erwise non-radiant polariton states via emission of bulk 
substrate photons. It was also studied how the radia- 
tive decay properties of quantum wire excitons get modi- 
fied when the wire is embedded in a microcavity, that is, 
via emission of cavity photons (Refs. dHH and references 
quoted there) . As one-dimensional arrays are considered 
for the directed energy transfer applications, their inter- 
action with the environment may also be used to achieve 
certain purposes such as, e.g., to counteract losses by 
embedding the array in the gain medium^ 

In this paper we discuss ID coherent dipole excita- 
tions that are formed in the neighborhood of the planar 
reflecting substrate in the range of frequencies to where 
the dielectric constant of the substrate medium e(u) < 
and the substrate does not support bulk photon modes. 
This is the situation that is most easily implementable 
in the vicinity of a metallic surface and which, in fact, 
has recently received a considerable attention in the con- 
text of both organic excitonsii and dipole excitations of 
nanoparticles (Ref. [l8| and references therein). The sub- 



strate in general affects both the dispersion of the exci- 
tons and their life-time. We will show that the presence 
of the substrate may result in new interference patterns 
and leads to a plethora of behaviors depending on the 
polarization and frequency of the excitation as well as 
on the distance from the interface. For the decay width, 
this is illustrated in Fig. Q] that exemplifies substantial 
differences with the vacuum result |T]) and is discussed 
in more detail later. This Figure demonstrates not only 
a modification of the exciton decay into vacuum photons 
but also the decay into substrate surface plasmon (SP) 
modes (that occur for e(oj) < —1), significance of the 
latter channel strongly increasing upon approach to the 
interface. 

Importance of the coupling of an individual dipole ex- 
citation to substrate SP modes and of associated reso- 
nance decay and scattering phenomena have already been 
stressed both for molecule o 19 ' 20 ' 21 and nanoparticles^ 
This coupling strongly increases as the exciton transition 
frequency approaches the "resonance" region of the sub- 
strate in which e(ui) is close to — 1. The increase is, of 
course, what should be expected from the theory of elec- 
trostatic image forces^ 2 - which features the combination 
factor 

Qo = (1 - e)/(l + e) (3) 
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FIG. 2: As in Fig. [T] but for d = 50 a and e = —1.2. Note that despite the separation distance here is 5 times larger than in 
panels (d-f) of Fig. [T] the decay into substrate surface plasmons is much stronger. 



for the magnitude of image charges. One should be 
aware, however, that this is also the region where both 
retardatio n 19 ' 20 and dissipation 1 - 1 - effects are particularly 
important. The corresponding enhancement of the de- 
cay of ID excitons into SPs is seen in the illustration 
of Fig. [U (e = —1-2) where, in comparison with Fig. [1] 
(e = —10), it clearly is a dominant decay channel; con- 
sequently the fluorescence efficiency is greatly reduced 3 
One also appreciates the fact that the resonant enhance- 
ment, as shown, occurs over the already "super-radiant" 
vacuum decay rate. 

Both Figs. [T] and [2] have been calculated with negligible 
exciton scattering and substrate losses (see Sec. IIVI for 
discussion of dissipation effects), hence they are reflective 
of the full conservation laws for our exciton-photon-SP 
system. So excitons with wave vectors |<7| > k cannot 
decay into vacuum photons, while emission of SPs can 
occur only for \q\ < n p (k) where 

Kp = Me/(e + i)] 1/2 (4) 

is the well-known (e.g., Ref. [Toh wave vector of the SP at 
an appropriate frequency. In such an idealized system, 
exciton-polaritons with larger wave vector magnitudes 
would be non-emissive. An example of the corresponding 
qualitative picture of the dispersion spectrum of exciton- 
polaritons in the chain is shown schematically in Fig. [3] 
The non-emissive branch 2 there can be thought of as 
of excitations representing a coherent mix of the exci- 
ton, photons and SPs of the same momentum projection 
along the chain, the relative weight contributions depend- 
ing on this momentum. (The branch splitting exhibited 
in Fig. [3] would not take place for y-polarized excitons as 
consistent with the SP polarization.) We cannot exclude 
that surface plasmon guiding by chains of nanoparticles 
found in recent numerical simulations 2 ^ is related to the 
formation of the discussed bound exciton-SP states. 

As the SPs are surface states characterized by 2D wave 
vectors, the inverse square-root singular behavior in the 
decay rates of ID excitons upon q — > k p seen in Figs. Q] 
and [2] has the "dimensionality" origin similar to the one 
taking place for 2D (quantum well) excitons decaying 
into 3D vacuum photons 1 ^ 24 ' 25 and ID (quantum wire) 
excitons decaying into 2D cavity photons. 8 
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FIG. 3: Schematically (not to scale!), an example of the pos- 
sible idealized dispersion of x-polarized ID exciton-polaritons 
in the form of their frequency u) = E(q)/h as a function of 
the wave vector q. Here the exciton-polariton spectrum (solid 
lines) is split into two branches: the states of branch 1 decay 
via emission of photons and/or SPs, the states of branch 2 
are non-emissive. This spectrum is a result of the interaction 
of the bare excitons with vacuum photons, whose spectrum 
is reflected in the limiting long-dash line, and with substrate 
surface plasmons, whose spectrum is reflected in the limiting 
short-dash line. ui sp is the asymptotic SP frequency. 



In what follows we will elaborate on the interaction 
of ID excitons with the neighboring reflecting substrate 
using the simplest model of a chain of molecules with a 
single molecular dipole transition (Frenkel-like excitons) . 
We note that some resulting features would be generically 
valid for Wannier-Mott excitons as well. We however 
do not pursue here an explicit analysis of the Wannier- 
Mott excitons that would require specific modifications 
for their bare dispersion as well as the influence of the 
substrate on the exciton binding and oscillator strength. 
While we chose, for certainty, to use an example of the 
metallic substrate and the corresponding surface plasmon 
excitations, it should be clear that the same effects would 
occur for a dielectric substrate whenever it can feature 
the region of negative e(uj) and the surface polariton ex- 
citations. 
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II. FRENKEL EXCITON-POLARITONS IN A 
MOLECULAR CHAIN. VACUUM RESULTS. 

Renormalization of a single-molecule transition of en- 
ergy Ef and the formation of an electric-dipole exciton 
band in a molecular aggregate can be derived in various 
frameworks (see, e.g., Refs. fll|26l) with identical results. 
Here we will use a simple and physically transparent de- 
scription having a clear underlying semi-classical anal- 
ogy. In the Heitler-London approximation (energy Ef is 
much larger than all other energies involved), the exci- 
ton Hamiltonian for interacting identical molecules can 
be written as 

H = {Ef + V ) B nBn + J2 V n - m {B\B m + h.C.) 
n n > m 

= Y,E{q)B\B q , 



where B\ creates an excitation on the nth molecule 
whereas W q = N^ 1 / 2 x ^ n e l i na £}t creates an excitation 
of the wave vector q in a chain of N molecules with the 
intermolecular spacing a. Correspondingly, the exciton 
energy 



electric fields; various aspects of it have been used for 
different geometries (e.g., Refs. [TollTllll^ and references 
therein) . In this paper we will not pursue solving specific 
self-consistent problems that may depend on a multitude 
of parameters and, instead, be discussing the self-energy 
for a given real value of parameter k (that is, the os- 
cillating frequency) for different values of the excitation 
wave vector q. For our numerical illustrations in this 
paper we chose a representative value of k correspond- 
ing to the wavelength A = 2n/k = 400 a, a reasonable 
magnitude for the optical region of the spectrum given 
typical spacing a ~ 10 A. Understandably, typical values 
for nanoparticle systems would be different X If neglect- 
ing the retardation effects (purely electrostatic fields), 
the value of k would have to be set equal to zero. 

Let us briefly review the application to a molecular 
chain in vacuum (see also Refs. and 0) ■ Consider the 
standard 22 retarded electric field at the point r = r r 
from the oscillating point dipole in vacuum: 

E v (fc,r) = 2^{fc 2 [p-f(f. p )] 

^-7) W-P)-P]} (9) 



E(q) = E f + Vq + 2^ V n cos(qna), (5) an d the corresponding 



n>0 

where Vq represents a possible renormalization for a sin- 
gle molecule 3 -^ while V n the effective intermolecular in- 
teraction as mediated by the electromagnetic field. The 
interaction is handily expressed via the semi-classical ex- 
pression 



Vn 



p ■ E(nax) 



(6) 



with p being the molecular transition dipole moment and 
E(nax) the electric field produced by the dipole p = pp 
at the distance na along the chain axis chosen to be along 
x (we use caps to denote unit vectors). 

If the electric field E is the total retarded field, how- 
ever, it is in fact also a function of the oscillating dipole 
frequency lo = ck and has both real and imaginary parts. 
Equation ((SJ) then has to be rewritten as a more general 
equation 



S v (fc,r) = -p-E v (fc,r). 

This expression turns out to be directly applicable even 
for the decay of a single molecule (see also Ref. Hoh : in- 
deed, calculating —2 Im{E v (fc, r — > 0)} immediately leads 
to the well-known decay width To derive the decay 
rate for a ID exciton state with wave vector q, one aug- 
ments this decay by the sum of contributions from other 
molecules: 

T v (fc, q) = r°(fc) - 4 cos(qna) Im{S v (fc, nax)} 

n>0 

resulting, after evaluation of the sum, in already quoted 
Eq. ||TJ). Exemplifying a general feature of self-energy 
corrections, a direct inspection easily verifies that, for a 
fixed k, Eqs. (p} and (J2J> satisfy, as expected, 



E(q) =E'- iT/2 = E f + E(ft, q), 



(7) 



J2r v (k,q)=NT a (k). 



(10) 



involving the self-energy correction S(fc, q), the real-space 
transform of which 



With the real part of the field ([9|), one immediately 
obtains the effective resonant interaction matrix element: 



S(fe, ar) = -p • E(fc, xx) 



(8) 



1 



serves to replace ©. Equation self-consistently (via 
k = E/Tic) determines both the dispersion E'(q) = Ef + 
Re{E(/c,o)} and the decay width T(q) = -2Im{S(fe, q)} 
of the exciton-polariton states as a function of their wave 
vector q. This simple approach, alternatively formulated 
in terms of Green's functions, is both physically appeal- 
ing and powerful as it involves only classically calculable 



Re{Y, v (k,x)} = (1 -3 cos 2 



- ( 1 - cos 2 



cos(fca;) ksin(kx) 



k 2 cos(kx) 



(11) 



exactly the same result that would be derived in the pic- 
ture of the virtual photon exchanged Calculating the 
corresponding sums with n > (and disregarding the 
irrelevant single-molecule renormalization 31 ) for many 
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molecules in the long wavelength, go < 1, expansion, 
one arrives at 



Re{E v (fc,g)} 



1-3 cos 2 



4C(3) 



3<? 2 



(fc 2 - q 2 ) b] + (1 - cos 2 6) k 2 b, (12) 



where C(3) ~ 1.202 and b = In (|fc 2 - g 2 |a 2 ). (It is useful 
to note that approximation (| 12(1 actually works very well 
over a sizable portion of the exciton Brillouin zone.) The 
known logarithmic divergence in Eq. (| 1 2[) upon q — > fc 
signifies the splitting of the exciton-polariton spectrum 
into two branches^" 7 - as caused by the radiative zone 
component of the field at cos 9 ^ 1 in Eq. (fTTj) . Of course, 
corresponding to this divergence there is a non- vanishing 
decay rate at q — > k in Eq. [TJ The vanishing of the latter 
takes place only at cos 6 = 1, and in this case the exciton 
dispersion exhibits a single-branch behavior with a cusp 
at g = fc. 

The electrostatic part of Eq. (112|) features a non- 
analytic behavior cx q 2 ln(ga) 2 at \q\ ^> k due to the 
long-range nature of the dipole-dipole interaction mak- 
ing the exciton dispersion "steeper", the behavior that 
recently attracted attention in the context of exchange- 
interaction effects for excitons in carbon nanotubes. 27 
The overall width of the bare exciton zone as seen in 
Eq. (JT2J) is scaled with the energetic parameter 



J = p 2 /a 3 



(13) 



establishing the unit for the nearest-neighbor electro- 
static dipole-dipole interaction. To appreciate the scale 
of energies involved: with p = 1 Debye and a = 10 A, 
for instance, J ~ 0.014 eV. With reasonable variations 
of values of p and a, J could reach magnitudes ~ 0.1 eV. 



III. MOLECULAR CHAIN NEAR THE 
INTERFACE. 

In the vicinity of the interface with a metallic/dielectric 
body, the total electric field can be conveniently repre- 
sented as a sum of the primary, vacuum, field, discussed 
in section [U and the secondary field due to the induced 
response of that body: E = E v + E s . Correspondingly, 
the self-energy of the exciton-polaritons is also repre- 
sented as E = E v + S s . In what follows we discuss the 
induced contribution E s (fc,g) coming from a half-space 
characterized by the dielectric constant e = e(u>) taken 
at the frequency uj — ck in the geometry of the molecu- 
lar chain (along x) being parallel to the separating inter- 
face (xy-plane) at the distance zq = d/2 from it (d is the 
distance between a dipole and its image). 

The problem of an electric dipole near a metal- 
lic/dielectric half-space is a classical problem first treated 
by Sommerfel d 10 ' 28 and whose solution is available in dif- 
ferent forms. Here we find it convenient to adopt the 
expressions for the electric field as derived in Ref. [29J. In 
the context of our application for fields along the chain, 



it matters how the dipoles are oriented with respect to 
both the chain and the interface. For a chain of electric 
dipoles of an arbitrary polarization p, one easily finds 
that 

E s (fc, a) = p 2 E^(fc, a) + p 2 E^(fc, x) + p 2 E|(fc, x), 

where axes-related components can be rewritten from re- 
sults in Ref. 29 as follows: 



^T, s x (k,x) = r Kdne-^ d f^J-( 
P Jo V 2 



for x-polarized dipoles, 



Kx) H J + (kx) 

2 7 

(14) 



1 f°° I ' iO k 2 P 

-jrEtfoa;) = / ndne-i d -^-J+Ux) + -—J-Ux) 
P Jo V 2 2 7 

(15) 

for y-polarized dipoles, and 

1 f°° K 2 (l 

^T,l(k,x)=J Kd K e-~< d — -Jo(kx) (16) 

for z-polarized dipoles. In Eqs. (ll4lll6|) . 

J±(x) = J (x) ± J 2 (x) 

are composed of Bessel functions of the first order while 
parameters 



Q = 



and 



7e ~ £7 



7 = (« a - fc 2 ) 1 / 2 , 7e = _ e e) 



p = 



7e +7 



2N1/2 



(17) 



(18) 



(for negative u < 0, u 1 ^ 2 — — i(— u) 1 / 2 should be used in 
Eq. I|18[l .) One can straightforwardly verify that the no- 
retardation limit (fc = 0) of the above expressions leads 
to usual electrostatic fields of image dipoles. 

Representations (114H16[) are quite meaningful physi- 
cally. One immediately observes that the pole in the 
parameter Q in (|17p occurs at k equal to k p in Eq. Q, 
the wave vector of the surface plasmon (surface polari- 
ton) whenever real e < — 1. Together with p-polarized 
photons, SPs make the Q-containing contributions in 
Eqs. (|14j|16p . The parameter P-containing terms, on 
the other hand, correspond to contributions to the elec- 
tric fields from s-polarized photons. Expressions (|14j|16|) 
taken at the source point x = would describe the effect 
of the half-space on the electronic transition in a single 
molecule as studied in Refs. Il9ll20l . 

We now turn to the g-dependent self-energy of an ex- 
citon in a long chain of molecules. Restricting to the 
chain-to-interface distances larger than the intcrmolec- 
ular spacing: d 3> a, one can safely use a continuum 
description of the half-space response: 



E s (fc,g) = 2 



dx 

— E s (k,x)cos(qx). (19) 
a 
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FIG. 4: Real (upper panels) and imaginary (lower panels) of the self-energy E(fe, q) of exciton-polaritons as functions of the 
reduced wave vector q/k for a fixed value of k = 27r/400a, e = —1.1 and d = 10a. Self-energy is shown with respect to J 
from Eq. (|13[) . The vertical short-dash lines show the position of the surface plasmon wave vector, K p /k. In the upper panels, 
the long-dash lines show the electrostatic (no-retardation) results, dash-dotted lines the electrostatic results for the exciton 
dispersion in vacuum. 



As with the vacuum case, all results are clearly even 
functions of q; to simplify expressions, we will therefore 
continue assuming q > 0. Transformation (|19[) for indi- 
vidual dipole contributions p4H16[) is facilitated by the 
integrals^ 

I dx cos(qx)J {nx) = — ^-rjrQ{n-q), (20) 

o (K - T) 1 



1 - 2q 2 In 2 

dx cos(qx),h(nx) = — -5 5-7 -jtt 6(k - q) (21) 

(«r — q ) * 



so that Eq. |[T6|) . e.g., is transformed into 



4X? B ( A)9 )=2 T-^ 
P A 7(«- 



K 3 Qe- 7d 



^2)1/2 



(22) 



and similarly for Eqs. (fT4")) and (jTSJ) . Step- functions in 
Eqs. (|2TJ|) and (|2T]l . as is also reflected in Eq. (|2"2"j). have a 
clear physical significance of the energy-and-momentum 
conservation limitation. 

It is also useful and meaningful to note the no- 
retardation limit (k = 0) of the above expressions, when 
Q in Eq. (fl7|) becomes equal to the electrostatic combi- 
nation ([3]), and the exciton dispersion would be affected 
by the image dipoles as 



poo -Kd 

E s x (q) = 2q 2 — 

hP J q (^-q 2 ) 1 ' 



— 2 dK = 2q 2 K (qd), 
(23) 



Qop 2 

and 



^ y (q) = 2 J°°{n 2 - q^e-^dK = ^(qd), 

(24) 



S s z (g) = E s x ( g )+S s y ( g ), 



(25) 



where Kq(x) and Ki(x) are the modified Bessel func- 
tions. 

The effects of these real image corrections are displayed 
in Fig. 2] showing an example of the total g-dependent 
exciton self-energy E(fc, q) including both vacuum and 
secondary field contributions. The figure is a result of 
an illustrative calculation for an idealized (no dissipa- 
tion) system in the resonance region (e = —1.1 so that 
k p ~ 3.3fc) at a relatively close (d = 10 a) distance from 
the interface - to demonstrate a possible magnitude of 
the effects. The figure shows both the retardation effects 
as well as the fact that upon the increase of the wave vec- 
tor q, the exciton dispersion approaches the electrostatic 
behavior. It is clear that, in principle, the electrostatic 
image forces may appreciably modify the exciton disper- 
sion. 

In an idealized system, the secondary field contribu- 
tions to the imaginary part of the self-energy (exciton 
decay) can come only from two sources. First, it is the 
regions of variable n in the integrals where the parameter 
7 (|T51) is imaginary. As Eq. (f2"2"]l shows, that happens only 
for q < k - this is the source of the modification by the 
substrate of the exciton decay into vacuum photons that 



7 




q/k q/k q/k 



FIG. 5: Decay width V(k, q) of the exciton-polariton as a function of the reduced wave vector q/k for a fixed value of k and 
e = —0.8. Two rows correspond to different dipole-to-image-dipole distances: upper panels (a-c) are for d — 600a, lower panels 
(d-f) for d = 10a. There is no surface plasmon channel here, and the total decay width is due to the vacuum photons only. 



wc explicitly illustrated in Fig. [TJ Second, and dominat- 
ing in Fig. |U source is the surface plasmon pole k = n p 
in the parameter Q in the integrals. As Eq. (|2"2")1 shows 
again, this pole contributes only for q < n p signifying the 
decay of the exciton into a SP. The pole contributions to 
— Im{S(fc, q)} are immediately calculable and can be rep- 
resented by a combination of a common factor 



47T 



|e|V2(i 



2 4 

p k p 

ak 2 



exp 



(26) 



and polarization- and q-dependent co-factors: 

q 2 /K p (n 2 p - q 2 ) 1/2 (27) 



for x polarization, 



(«2 - « a ) 1/a /»p 



for y and 



|eK/(4 - S 2 ) 1/2 



(28) 



(29) 



for z. 

Accompanying decay's inverse square-root singularity 
in Eqs. ([STJ) and (|29"|). also seen in Fig.[4|d) and (f), are 
the diverging discontinuities in the real parts (panels (a) 
and (c) of Fig. Ql of the type familiar from the studies of 
the decay of 2D excitons into 3D vacuum photons i 1 ' 5 ' 24 ! 25 
Due to such a discontinuity, the resulting self-consistent 
dispersion of, e.g., x-polarized excitations would split 
into two branches as shown in Fig. [3J On the contrary, 



the decay rate of y-polarized excitons vanishes at q — > k p 
in Eq. (|2"8)) . This is a consequence of the polarization 
of SPs whose electric field can have only longitudinal 
and perpendicular to the interface components (like p- 
polarized photons). As q — > k p , the surface plasmons 
would be emitted along the chain direction, hence their 
field would have no y-components to interact with y- 
excitons. (There is a similarity here with the decay of x- 
polarized excitons in vacuum, whose decay rate vanishes 
at q — > k, Eq. JJJ.) Correspondingly, the real part of the 
self-energy in Fig. [31(b) does not exhibit a discontinuity 



at q 



(the discontinuity is in the derivative) with 



the resulting exciton spectrum consisting of one branch 
only. Quite clear is also cx q 2 behavior for x-excitons in 
Eq. ([2T| . it has the same origin as in the electrostatic 
effect ([23| - at q — ► 0, the exciton polarization of the 
chain becomes uniform and there would be no polariza- 
tion charges (a vanishing spatial derivative) to induce 
images in the substrate. 

The SP excitations in the substrate exist in the region 
of frequencies where e(o>) < — 1. Qualitatively different, 
therefore, for a reflecting substrate is another region of 
frequencies in which —1 < e{uS) < 0. In this case exciton- 
polaritons in the molecular chain can decay only into 
vacuum photons. The (idealized) substrate then serves 
as to modify the dispersion of excitons and their radia- 
tive decay. An example of the radiative decay modifi- 
cation is shown in Fig. [5] calculated for e = —0.8 and 
to be compared with Fig. [TJ Both figures feature the 
same distances from the interface but different magni- 



8 



tildes and signs of the image charges, Eq. ([3]). These 
two facts explain similarities and differences in the ra- 
diative decay patterns for the two examples. Among the 
important common features, one should notice "oscillat- 
ing" g-dependences of the decay rate at larger distances 
and the vanishing of the radiative decay at q — > k. As 
the chain is moved further away from the interface, even 
more undulations would be observed in the g-dependence 
of the decay rate as a result of the interference with the 
image dipoles. 

IV. EFFECTS OF DISSIPATION IN THE 
SUBSTRATE 

Real metallic (dielectric) substrates are characterized 
by finite dissipation (losses) that are taken into account 
phenomenologically via the imaginary part of the dielec- 
tric function: 

e(w) = e'(co) + ie"(w). 

In the context of the effects of e" on SPs, for a given 
real frequency u), wave vectors of the SPs consequently 
acquire imaginary parts as welli^ 

ftp — ftp + iftp, 

where ftp describes the damping of the SP modes. An- 
other way to look at of the uncertainty (broad- 
ening) of the SP's wave vector. Hence the momentum 
conservation law in the exciton decay into a SP would 
not be obeyed exactly as in idealized systems discussed so 
far. One thus should expect a broadening and extension 
of the range of finite exciton decay rates into the region 
of the exciton wave vectors q > ft p . Moreover, even with 
relatively small losses, e" <C |e'|, the relative damping of 
the SP can become substantial in the resonance region of 
e' close to — 1. Indeed, in this case^ 

<K - 4tY) ' (30) 

and the denominator in the r. h. s. of Eq. ([50)) could 
"compensate" for the smallness of e" and completely de- 
stroy the notion of a coherent SP. 

Figure [6] shows results of the calculation of the exci- 
ton decay for the substrate with the dielectric constant 
e = —1.2 + 0.05« and variable distances to the inter- 
face. As compared to the idealized system, one immedi- 
ately notices very broad, especially at smaller distances, 
q 3> ftp regions of the exciton decay. For this particu- 
lar value of the complex dielectric function, however, the 
"uncertainty" ratio (|5D|) is only about 0.1. Thus the ex- 
tremely broadened g-region of the exciton decay is not 
due to the broadening of the SPs themselves, the latter 
can explain only the formation of the finite magnitude SP 
peaks clearly seen in Fig. [5J Of course, with complex di- 
electric constants, the imaginary part of the exciton self- 
energy discussed in Sec. IIIH is contributed to, in principle, 



by the whole integration range in integrals like Eq. (|22[) . 
rather than just by the around-the-pole region. What is 
reflected in the broad "wings" in Fig.[S]is the result of the 
ordinary "incoherent" Joule losses due to the oscillating 
electric field of the exciton polarization penetrating into 
the substrate. Especially illuminating in this regard is a 
second maximum on the top- most curve in Fig.[5Ja). In- 
deed, similarly to the already discussed in Sec. lIIII image- 
charge effects, the electrostatic component of the electric 
field of x-polarized excitons has a q- and d-dependence 
reflected in the r. h. s. of Eq. (|2"5|) , the maximum of which 
is achieved at q ~ 1/d. As the distance from the inter- 
face increases, the role of such short-range energy transfer 
from the higher-g excitons to the substrate decreases and 
one can "tune" it off while still having an appreciable de- 
cay rate of the lower-g states into SPs. In addition to the 
largest-distance curves in Fig. [6l this point is illustrated 
in Fig. [7] showing clear SP emission effects both close to 
and away from the resonance region. The distance de- 
pendence of the SP emission intensity is governed by the 
factor f2^1) and one could likely "optimize" the relation- 
ship between the coherent and incoherent energy transfer 
based on the system regime parameters. 



V. DISCUSSION. 

The modification of the electric fields in the presence 
of the reflecting substrate leads to substantial changes 
of the properties of ID exciton-polaritons in a dipole- 
coupled chain above the substrate. The substrate mod- 
ifies the interaction of the exciton with vacuum photons 
and, in addition, can engage a new and efficient interac- 
tion channel - with surface excitations of the substrate. 
While these statements also apply to a well-studie d 19 i 20 
case of a single molecule near the interface, the cooper- 
ative interaction of many molecules in the chain and of 
their images may lead to qualitatively and quantitatively 
interesting effects. In this paper, we discussed such ef- 
fects for coherent delocalized ID exciton-polaritons that 
are well characterized by their wave vectors q. We em- 
phasize that it is actually individual g-states that would 
exhibit qualitatively different coherent properties - if av- 
eraged uniformly over all q, a (near) reduction would 
occur to the effects of the substrate on an individual 
molecule (similarly to discussed Eq. (|10|) for vacuum). 
Immediately noticeable in our illustrative examples, for 
instance, is a difference of the behavior of x- and y- 
polarized ID exciton-polaritons - while those cases of 
the molecular transition dipoles parallel to the interface 
would be equivalent for single molecules. Properties of 
ID exciton-polaritons depend specifically on their polar- 
ization, frequency with respect to the substrate dielectric 
dispersion and distance from the interface. We reiterate 
that, while our illustrations used parameters appropriate 
for molecular systems, qualitatively similar effects are ex- 
pected for other systems such as chains of nanoparticles. 

A common feature for all exciton polarizations and 
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FIG. 6: Decay width V(k, q) of exciton-polaritons as a function of the reduced wave vector q/k for a fixed value of k — 27r/400 a 
in units of vacuum F v = 2irp 2 k 2 /a and for different (indicated in the panels) exciton polarizations. Here substrate's dielectric 
constant is complex: e = —1.2 + 0.05i. Each of the panels contain 4 curves corresponding to different distances between the 
molecular chain and the substrate (top to bottom): d = 10 a, 20 a, 30 a and 50a. 
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FIG. 7: As in Fig.[6]but for d = 50a and two different dielectric functions (top to bottom): e = — 2.0 + 0. li and e — — 10.0 + 0.5i. 



conditions considered (see Figs. QJ El E]) is that the rate of 
the radiative decay into vacuum photons in the presence 
of the substrate vanishes at q = k = co/c, that is, where 
the exciton dispersion curve would cross the vacuum pho- 
ton dispersion line (see, e.g., Fig. [3]) - this is distinctly 
different from the case of the chain in the vacuum where 
such vanishing would occur only for excitons polarized 
parallel to the chain. The consequence of this effect of 
the image dipoles is that no branch splitting occurs at 
q = k. At distances from the interface comparable to 
the wavelength 27r/fc of vacuum photons, the interfer- 
ence with the radiative fields of the image dipoles results 
in undulated patterns of the radiative decay (panels (a-c) 
of Figs. [TJ [5]) as functions of q; the larger the distance the 
more undulations would take place. 

It is important to note that, as a result of integration 
(flU]) over many molecules in the chain, the magnitudes of 
these and some other g-dependent features fall off with 
the distance from the interface slower than their aver- 
ages characteristic of single molecules. Consider, e.g., 
the electrostatic exciton energy shifts due the interaction 
of the molecular dipole(s) with their image(s). For a sin- 
gle molecule, the corresponding transition energy shift is 
~ p 2 /d 3 . The maximum of g-dependent shifts, however, 
as Eqs. (|23M25I) show, would be ~ p 2 /ad 2 - that is, much 
larger at distances d>o (not in the immediate proxim- 
ity of the interface) . This may provide an opportunity of 



easier experimental identification of such shifts than for 
individual molecules. 

While the distance dependence (|2"o| of the decay rate 
of the excitons into surface plasmons is the same as for 
single molecules, g-dependent enhancing factors (|27|l and 
(|29p. seen as the peaks in our illustrations, may facilitate 
better experimental verifications. As the rate of the emis- 
sion of SPs is g-dependent, we speculate that the chain 
excitons could perhaps serve as directional sources of SPs 
- in expectation of a more efficient emission along the 
chain for x- and z-polarized excitons and perpendicular 
to the chain for y-polarized excitons. We recall that these 
polarization assignments stem from the fact that the elec- 
tric field of the SPs lies in the plane made by the wave 
vector and the normal to the interface. As the small-g 
exciton decay into SPs can be strongly enhanced in com- 
parison with the decay into vacuum photons, it is likely 
that the inverse process of the excitation of the chain 
by SPs could also be exploited. Non-emissive exciton- 
polaritons with larger q may also present an opportunity 
to be used for SP guiding. 23 

Various scattering and dissipation processes are known 
to be able to strongly affect features characteristic of ide- 
alized systems. We particularly discussed the chain exci- 
ton quenching due to the "incoherent" energy transfer to 
the substrate - and this does not exhaust the list. It suf- 
fices to also mention, e.g., the scattering by phonons in 
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molecular chains or Joule losses in the chains of metallic 
nanoparticles. The fast scattering between different in- 
states would result in the thermalized population of the 
exciton-polaritons so that observed decay rates are ther- 
mally averag ed (see Refs. 0, i, M 

for some specific ID 
applications). We therefore presume that the best condi- 
tions to experimentally address finer g-dependent effects 
we discussed in this paper would be low-temperature 
spectroscopic measurements. 
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